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has  been  su.'igostod.  The  ^00  and  GOO  l:m  discontinuity  in  the  upper  mantle  have 
been  rorrelated  with  the  shir()  chatujes  in  travel-time  gradients  at  1900  and 
2890  I;m,  respectively. 

The  second  problem  deals  with  the  morphology  of  the  Benioff  zone  in 
southcenti'al  Al(;ska.  The  results  show  tfiat  the  underthrust  North  Pacific 
lithospheric  plate  extends  north  up  to  about  G3^N  across  the  Alaska  Range.  Also 
the  plate  shows  a possible  break  under  the  Yentna  j'iver  Prince  Wiliam  Sound  area, 

"^he  third  problem  deals  with  the  P-wave  first-motion  solutions  of  small 
magnitude  earthquakes  located  in  central  Alaska.  The  earthquakes  with  focal 
depths  less  than  GO  km  show  either  strike-slip  or  normal  faulting  focal 
mechanism,  while  those  with  focal  depths  greater  than  the  above  value  show 
nearly  pure  dip-slip  motions.  The  solutions,  in  general,  show  a northwest- 
southeast  trend  of  the  pressure  axes.  Also  the  regional  features  of  the 
focal  mcchanif'’,  show  good  correlation  with  those  pi'edicted  from  North  Pacific 
plate  tectonics. 

The  fourth  problem  deals  with  the  surface  and  body  wave  magnitude  relations 
for  earthquakes  and  explosions  for  tv.'o  regions,  Nevada-California  and  central 
Aleutian.  In  the  [■‘,s-m.  relations,  the  earthquakes  from  both  the  areas  show 
clear  separation  from  the  explosions. 
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1.  TRAVEL-TIME  RELATIONS  TOR  THE  CRUST  AND  UPPER  MANTLE  P-WAVE  PHASES 

FROM  CENTRAL  ALASKAN  DATA 

By  N.  N.  Biswas  and  B.  Bhattacharya 
ABSTRACT 

This  section  of  the  report  deals  with  the  travel -time  relations  for 
the  crust  and  upper  mantle  P-wave  phases  in  the  epicentral  distance 
range  of  100  to  4500  km.  The  short-period  data  analyzed  were  from 
underground  nuclear  explosions  in  the  central  Aleutian  and  Nevada  regions, 
from  earthquakes  located  in  central  and  southcentral  Alaska  and  along 
the  nortfi  Pacific  seismic  beltsi  the  data  were  recorded  by  the  interior 
Alaskan  seismic  network. 

In  the  epicentral  distance  range  of  100  to  600  km,  the  first  arrival 
P-wave  phases  are  P*  and  P^;  the  former  has  been  identified  as  the 
earliest  arrival  for  a<300  km  with  an  apparent  velocity  of  7.05  km/sec 
while  for  epicentral  distances  greater  than  300  km,  the  latter  phase 
arrives  first  and  slows  an  apparent  velocity  of  8.07  km/sec.  A tenta- 
tive Mohe  depth  of  4C  km  for  central  Alaska  has  been  obtained  from  this 
study.  This  value  is  in  close  agreement  with  that  derived  for  southern 
Alaska  by  others. 

The  data  for  the  epicentral  distance  range  of  700  to  4500  km  could 
be  explained  by  three  linear  relations:  the  first  extends  from  700  to 
1900  km  with  apparent  velocity  of  8.29  km/sec,  the  second  lies  between 
1900  and  2850  km  with  apparent  velocity  of  10.39  km/sec,  and  the  third 
extends  from  2850  to  4500  km  with  apparent  velocity  of  12.58  km/sec. 

In  the  epicentral  distance  range  of  1900  to  2600  km,  the  first  P-wave 
ai'"ivals  from  earthquakes  in  the  Aleutian  area  were  found  to  be  systematically 
late  by  3 to  4 sec  compared  witn  those  from  central  Aleutian  underground 
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nuclear  explosions.  This  discrepancy  in  travel  tines  has  been  interpreted 
as  due  to  mislocatiotis  and  errors  in  origin  time  for  earthquakes  in  the 
above  area.  Also,  the  well  known  400  and  600  km  discontinuities  in  the 
upper  mantle  have  been  correlated  with  the  sharp  changes  in  travel-time 
gradients  at  1900  and  2850  km,  respectively. 

INTRODUCTION 

The  object  of  this  study  was  to  derive  the  regionalized  travel-time 
relations  for  the  crust  and  upper  mantle  P-wave  phases  from  the  data 
gathered  by  the  local  seismological  network.  Because  of  the  unavoidable 
source-receiver  disposition,  most  of  the  travel  paths  considered  for  the 
upper  mantle  phases,  except  for  P^,  traverse  the  complex  northern 
portion  of  the  Circum-Pacif ic  seismic  belt.  For  the  basal  crustal  phase 
P*  and  upper  mantle  phase  P^,  the  earthquakes  selected  lie  within  the 
perimeter  of  the  local  seismological  network.  Due  to  large  scatter  in 
travel-time  data  for  short  epicentral  distances  (a<100  km)  we  could  not 
study  the  intermediate  crustal  phase  P^.  We  anticipate  resolving  the 
characteristics  of  this  phase  from  chemical  explosions  planned  to  be 
carried  out  by  the  United  States  Geological  Survey  sometime  next  year. 

The  events  considered  for  this  study  which  lie  along  the  Alaskan 
coastal  belt  on  the  south  and  southv/est  of  the  seismic  network  are 
located  along  the  well  known  deep  sea  trench- island  arc  system.  The 
predominant  feature  of  this  system  is  a lithospheric  plate  sinking  into 
the  upper  mantle,  which  is  characterized  by  a low-velocity  zone,  phase 
change,  and  rapid  velocity  variation.  The  seismicity  indicates  a variable 
dip  of  the  sinking  lithosphere,  namely  40°  to  50°  in  southern  Alaska- 
central  Aleutian  region,  50°  to  60°  on  the  west  of  Rat  Island  and  90°  at 
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the  western  end  of  the  Aleutian  arc  (Jacob,  1972;  Davies  and  Julian, 

1972;  Engdahl,  1973).  Farther  west  along  the  Kurile- Kamchatka  seismic 
zone,  the  underthrust  plate  strikes  in  the  northeast-southwest  direction 
with  a dip  of  the  order  of  50°  to  60°  due  northwest  (Fedotov,  1968). 

The  above  regional  geometry  of  the  plate  is  in  accord  with  that  derived 
from  the  relative  motion  of  the  North  Pacific  plate  with  respect  to  the 
American  plate  along  the  southern  Alaska-Aleutian  region  and  the  former 
and  the  Eurasian  plate  along  the  Kurile-Kamchatka  region  (McKenzie  and 
Parker,  1967;  Le  Pichon,  <968;  Stauder,  1968a,  b,  1972). 

Epicenters  of  earthquakes  with  focal  depths  between  100  and  200  km 
located  along  the  trench-arc  syste..  show  positive  correlation  with  the 
positions  of  active  andesite  volcanoes.  Le  Pichon  et  al.  (1973)  summarized 
important  results  concerning  the  interrelationships  of  various  features 
observed  along  ihe  plate  bounda>--!es.  Moreover.  se\  ?ral  autnors  (Sorrells 
et  al.,  1971;  Toksoz  et  al.,  1971;  Jacob,  1972;  Davies  and  Julian,  1972: 
Abe,  1972;  Sleep,  1973)  derived  models  for  the  sinking  lithosphere.  For 
the  central  Aleutian  arcs,  Jacob  (1972)  proposed  a model  in  which  the 
sinking  plate  reaches  a depth  of  250  km  from  the  surface  with  7 to  10 
per  cent  higher  P-wave  velocity  than  the  surrounding  medium;  the  thickness 
of  the  plate  was  shown  to  be  80  km.  Quite  similar  models  have  been 
considered  by  Davies  and  Julian  (1972)  ana  Sleep  (1973)  for  their  studies 
of  the  explosion  data  from  Amchitka  Island.  The  close  similarity  of 
seismicity  between  the  central  Aleutian  and  Kurile-Kamchatka  seismic 
zones  perhaps  indicates  the  structure  of  the  dipping  plate  under  the 

latter  area  as  similar  to  that  proposed  for  the  former  area  by  Jacob  and 
others. 
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The  source  region  on  the  southeast  of  the  seismic  networks  belong 
to  the  well  known  strike-slip  fault  systems  marking  the  contact  planes 
between  the  Pacific  and  American  plates.  Besides  the  Nevada  underground 
nuclear  explosions,  these  events  were  chosen  to  compare  the  P-wave 
travel  times  observed  in  central  Alaska  for  propagation  paths  outside 
the  deep  sea  trench-island  arc  system.  The  results  have  been  discussed 
in  later  sections. 

The  receiver  region  is  also  complex;  the  major  structural  trends 
south  of  64°N  as  shown  in  Figure  1,  are  predominantly  strike-slip 
faults  which  parallel  the  broad  arcuate  structure  of  the  Alaska  Range. 
Previous  geophysical  studies  (Woollard  cu  al.,  1960;  Hanson  et  al., 

1968)  indicate  the  presence  of  a mountain  root  under  the  above  range. 

The  summary  of  Alaskan  seismicity  (Terr.  1S70)  shows  the  extension  of 
the  Aleutian  seismic  belt  in  interior  Alaska  where  it  widens  and  spreads 
over  a region  of  about  750  km  across.  The  hypocenters  of  earthquakes 
decrease  in  depth  from  about  150  km  in  the  Cook  Inlet  area  to  less  than 
50  km  along  the  northern  and  eastern  margins  of  the  inland  active  seismic 
zone. 

Davies  and  Berg  (1973)  have  provided  evidence  of  a sharp  westward 
dip  of  the  M-discontinui ty  just  north  of  Cock  Inlet.  From  hypocenter 
distribution  of  earthquakes  in  central  Alaska,  Van  Wormer  et  al.  (1974) 
showed  the  termination  of  the  Benioff  zone  in  interior  Alaska  some  75  km 
north  of  the  Denali  fault.  They  also  found  the  possibility  of  a break 
in  the  subducting  lithospheric  plate  in  the  Yenta  River-Prince  William 
Sound  area.  Lahr  et  al.  (1974)  and  Lahr  (1974)  made  similar  studies 
with  improved  seismic  control  around  Cook  Inlet.  They  obtained  a westward 
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dipping  Benioff  zone  reaching  a depth  of  150  km  from  the  surface  in  the 
above  area.  The  seismic  zone  along  the  sinking  lithosphere  in  the  depth 
range  of  50  to  150  km  was  shown  to  have  a thickness  of  about  15  to  20 
km.  It  seems  that  further  reduction  in  location  errors  may  show  the 
true  thickness  of  the  seismic  zone  smaller  than  these  values  (Engdahl, 
1973). 

Carder  et  al.  (1967)  reported  the  seismic  body-wave  arrivals  in  the 
distanc.'  range  of  0°  to  27°  from  the  nuclear  explosion  LONGSHOT  on 
Amchitka  Island.  During  two  latexplosions  (MILROW  and  CANNIKIN),  a 
number  of  permanent  and  semi-permanent  seismic  short-period  stations 
were  operated  by  a number  of  federal  and  state  agencies.  Biswas  (1973) 
reported  the  distribution  of  the  observed  P-wave  travel-time  residuals 
in  the  Alaska-Aleutian  region  by  combining  data  from  all  three  nuclear 
explosions.  The  results  are  shown  in  Figure  2;  the  values  in  parentheses 
are  mean  residuals  (5t  = <t^  - tj^>),  where  t^  is  the  depth-el lipticity 
corrected  observed  travel  time  and  t^  is  the  corresponding  value  of 
Herrin  et  al.  (1968)  for  a given  epicentral  distance. 

The  time  residuals  indicate  that  for  propagation  from  the  central 
Aleutian  area,  the  P arrivals  are  systematically  early  in  central  Alaska 
to  the  east  of  approximately  151°W.  These  may  result  from  structural 
complexities  either  at  the  source  region  or  at  the  receiver  region  or  a 
combination  of  both.  Thus,  prior  to  any  further  analysis,  it  is  necessary 
to  separate  the  observed  residuals  into  major  components.  Studies  to 
determine  station  re'iduals  from  teleseismic  events  located  in  different 
azimuths  with  respect  to  central  Alaska  are  in  progress.  In  the  following 
sections,  we  have  discussed  the  first  arrival  P-wave  travel-time  relations 
in  the  epicentral  distance  range  of  100  to  4500  km  from  short-period 
data  obtained  by  the  central  Alaskan  seismic  network. 
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RESULTS 


A.  and  Crustal  P-wave  in  Central  Alaska 

Because  of  the  lack  of  explosion  data  in  large  parts  of  central 
Alaska  for  the  epicentral  distance  range  in  which  the  crustal  phases 
and  P^  are  observed  as  first  arrivals,  vie  have  based  our  study  to  date 
solely  on  impulsive  first  arrivals,  we  have  based  our  study  to  date 
solely  on  impulsive  first  arrival  P-wave  phases  from  earthquakes  (1967 
through  1972)  distributed  inland  and  on  coastal  sections  of  central  Alaska. 
The  locations  of  the  epicenters  are  shown  in  Figure  3a;  the  results  have 
been  summarized  in  Figure  3b  and  3c. 

The  location  and  the  origin  time  of  the  earthquakes  considered  for 
this  study,  as  well  as  for  other  studies  discussed  in  later  sections 
were  taken  from  NOAA's  Monthly  Bulletins  (PDE).  The  observed  travel 
times  have  been  reduced  to  sea  level  by  applying  time  corrections  for 
focal  depth  and  receiver  elevation.  The  corrections  were  obtained 
either  directly  from  the  Herrin  table  (1968)  when  focal  depths  coincided 
with  the  table  values,  or  values  were  interpolated  from  the  above  table 
for  appropriate  focal  depths.  Finally,  the  ellipticity  correctior.t  were 
applied  to  the  depth-corrected  travel  times  from  Bullen's  tables  (1937). 

The  travel -time  data  show  a break  in  the  gradient  around  A = 300 
km.  This  can  be  seen  clearly  in  the  reduced  time  plot  (Figure  3c).  We 
thus  fitted  the  data  by  two  linear  relations  in  the  least-squares  sense. 

For  A<300km,  the  linear  relation  yielded  an  apparent  velocity  of  7.05 
km/sec  and  an  intercept  time  of  3.37  sec.  The  apparent  velocity  obtained 
for  the  P^  phase  (a<300  km)  is  8.07  km/sec  and  an  intercept  time  of  8.62 
sec.  The  first  P-phase  appears  to  represent  P*;  we  have  yet  to  establish 
the  values  for  P^  for  this  area  as  mentioned  in  the  previous  section. 
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Hales  and  Asada  (1966)  studied  P-wave  travel-time  relations  in  the 
epicentral  distance  range  of  80  to  400  km  in  southern  Alaska.  The 
sources  were  a series  of  chemical  explosions  located  in  the  Prince 
William  Sound  area.  Their  travel -time  relations  for  the  lower  crust 
show  an  apparent  velocity  of  6.98  km/sec  and  an  intercept  time  of  4.2 
sec.  These  agree,  within  the  experimental  accuracy,  with  our  results 
obtained  for  a<300  km.  Also,  from  the  analysis  of  data  obtained  along  a 
profile  close  to  the  Alaska-Canada  border.  Hales  and  Asada  (1966)  obtained 
an  apparent  velocity  of  8.07  km/sec  and  an  intercept  time  of  9.1  sec  for 
the  P^  phase.  Our  results  for  300  km  < a < 600  km,  though  obtained  for 
central  and  southcentral  Alaska,  are  very  similar. 

B.  Travel  Time  Relations  for  the  Upper  Mantle  P-Wave  Phases 

For  this  study,  short  period  oata  collected  from  1970  through  1972 
by  the  central  Alaskan  stations  have  been  analyzed.  Records  from  146 
out  of  about  340  earthquakes  were  found  with  clear  impulsive  first 
arrivals  and  the  arrival  times  could  be  read  within  + 0.2  sec.  The 
location  of  the  selected  events  are  shown  in  Figure  4,  and  regionwise, 
the  ranges  of  magnitude  and  focal  depth  are  given  in  Table  1. 

Earthquakes  with  magnitudes  less  than  the  lower  bound  of  Table  1, 
and  originating  on  the  western  side  of  Alaska  Peninsula  and  further 
west  along  the  island  arc  structure,  shew  a few  cycles  of  small  amplitude 
precursors  as  the  first  arrivals  followec  by  a large  amplitude  envelope. 
This  phenomenon  has  been  reported  before  by  Jacob  (1972)  for  the  LONGSHOT 
event  recorded  at  Point  Barrow.  Similar  features  were  observed  in 
western  Canadian  records  from  the  above  event  by  Davies  and  Julian 
(1972)  and  in  LASA  records  by  Davies  and  Frazer  (1970)  from  Aleutian 
earthquakes. 
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Jacob  (1972),  by  employing  ray- tracing,  suggested  that  the  large 
amplitude  envelope  could  be  either  a P-phase  refracted  from  the  top  of 
the  dipping  lithospheric  plate  along  the  Aleutian  seismic  zone,  or  the 
P-phase  refracted  from  the  600-km  discontinuity.  He  related  the  highly 
attenuated  first  arrivals  to  the  normal  P-phase  that  p-^opagates  just 
above  the  plate.  Davies  and  Julian  (1972)  made  similar  suggestions  for 
the  signal  characteristics  observed  in  western  Canada,  again  based  rn  a 
ray-tracing  method  for  the  LONGSHOT  event. 

We  have  examined  records  for  a large  number  of  Aleutian  earthquakes 
obtained  by  central  Alaskan  stations  having  high-  and  low-gain  recording 
channels.  Frequently,  the  high-gain  trace  shows  a distinct  impulsive 
first  break,  thereby  minimizing  the  reading  error  to  + 0.2  sec  while  the 
low-gain  trace  shows  a precursor  to  the  main  signal.  Our  results  discussed 
below  arc  based  solely  on  impulsive  first  arrivals. 

The  depth-el lipti city  corrected  travel  times  are  shown  in  Figure  5 
as  a function  of  the  epicentral  distance  in  the  range  of  700  to  4500  km; 
the  plot  of  the  reduced  travel  time  is  shown  in  Figure  6.  In  these 
figures,  we  have  also  included  travel-time  data  for  impulsive  first  P- 
arrivals  from  three  nuclear  tests  on  Amchitka  Island,  from  five  nuclear 
tests  in  Nevada  and  from  a f.^w  selected  earthquakes  located  in  the 
southeastern  azimuth  with  respect  to  the  central  Alaskan  seismic  network. 

It  should  be  noted  that  the  locations  and  origin  times  of  earthquakes 
were  obtained  using  the  Jeffreys-Bul len  tables,  whereas  those  values  are 
known  for  the  explosions.  Thus  the  travel  times  for  earthquakes  are 
subject  to  additional  uncertainties  because  of  errors  in  locations  and 
origin  times.  The  travel-time  values  corresponding  to  the  explosions 
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and  southeastern  events  are  shown  in  Figures  5 and  6 by  crosses  and 
hollow  triangles,  respectively.  We  incorporated  data  from  southeastern 
events  for  reasons  mentioned  in  previous  sections. 

A number  of  interesting  features,  as  discussed  below,  can  be  observed 
particularly  in  the  reduced  travel-time  plot  (Figure  6). 

As  the  epicentral  distance  increases,  the  scatter  in  the  data 
decreases.  This  feature  becomes  apparent  for  A > 2500  km  and  may  be 
explained  by  the  fact  that  the  errors  contributed  by  location  uncertainties 
or  the  fluctuations  in  travel-time  values  caused  by  structural  heterogeneities 
or  a combination  of  both,  become  small  compared  to  the  true  travel  times 
for  large  epicentral  distances. 

In  the  epicentral  dis.ance  range  of  2100  < a < 2500  km,  the  travel 
times  for  earthquakes  are  systematically  late  by  3 to  4 sec  compared  to 
those  from  explosions.  The  locations  of  the  events  with  respect  tc  the 
central  Alaskan  seismic  network  of  stations  for  the  above  distance  range 
extend  from  about  162“W  to  the  western  end  of  the  Aleutian  seismic  zone. 

It  is  well  known  that  the  location  of  LONGSHOT,  using  the  usual  location 
methods  with  teleseismic  data,  revealed  considerable  errors  in  computed 
values  of  the  focal  depth,  position  of  the  epicenter  and  origin  time. 

The  details  can  be  found  in  Tables  1 and  2 of  Lambert  et  al.  (Vol.  2, 

1970).  When  the  focal  depth  was  constrained  to  zero,  the  authors  found 
a northward  shift  by  20  km  in  the  computed  position  of  the  epicenter 
from  the  true  site.  Engdahl  (1973)  relocated  a number  of  intermediate 
depth  earthquakes  in  the  central  Aleutian  area  from  data  obtained  by  a 
local  seismic  network.  He  found  the  hypocenters  were  displaced  by  30  km 
upward  compared  with  those  determined  from  the  usual  location  method. 

An  error  in  epicenter  location  by  + 20  km,  or  in  focal  depth  by  + 30  km. 
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introduces  an  uncertainty  of  + 2 to  + 3 sec  in  travel  time  at  a 2200 
km,  the  epicentral  distance  conmon  for  the  Aleutian  events  from  central 
Alaskan  stations. 

In  view  of  the  above  results,  it  is  suggested  that  the  observed 
discrepancy  in  travel  times  between  earthquakes  and  explosions  in  the 
epicentral  distance  range  of  1900  to  2600  km  is  mainly  due  to  systematic 
errors  in  the  reported  values  (PDE:  location  and  origin  time)  for  the 
Aleutian  earthquakes.  This  interpretation  seems  necessary  but  may  not 
be  sufficient,  as  refraction  of  the  rays  emerging  in  central  Alaska  by 
the  sub-Aleutian  lithospheric  slab  may  partially  or  totally  invalidate 
the  great-circle  travel  path  assumption,  a factor  inherent  in  our 
analysis. 

In  the  epicentral  distance  range  of  700  to  1900  km,  19  reliable 
data  points  from  underground  nuclear  explosions  on  Amchitka  Island  were 
available.  These  data  have  been  supplemented  by  high-quality  data  from 
earthquakes  located  in  the  southeast,  south  and  southwest  directions 
with  respect  to  central  Alaska.  The  results  show  (Figures  5 and  6) 
that,  at  least  to  the  present  accuracy  of  our  measurements,  the  first- 
arrival  travel  times  for  the  above  epicentral  distance  range  cannot  be 
grouped  on  the  basis  of  the  direction  of  propagation  alone,  as  the  data, 
in  addition  to  following  a single  linear  trend,  exhibit  comparable 
scatter  from  all  nhree  azimuths.  However,  the  range  of  scatter  may  be 
reduced  by  incorporating  station  corrections  currently  not  available  for 
the  stations  used  in  this  study. 

In  the  epicenter  distance  range  of  about  3400  to  4000  km,  the 
travel-time  values  for  explosions  are  from  underground  nuclear  tests  in 
Nevada.  These  are  in  agreement  with  those  from  earthquakes  approximately 
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in  the  same  azimuth.  The  rest  of  the  values  between  2600  and  4500  km 
were  obtained  from  earthquakes  located  in  western  Aleutian  and  Kurile- 
Kamchatka  seismic  zones.  These  values  are  also  in  close  agreement  with 
those  from  the  southeastern  azimuth,  although  the  travel  paths  traverse 
dissimilar  tectonic  settings. 

The  travel-time  plot  as  shown  in  Figure  5 shows  three  linear  segments 
marked  by  two  sharp  changes  in  gradient,  one  around  2000  km  and  the 
other  arouno  2600  km.  These  can  be  seen  clearly  in  the  reduced  time 
plot  (Figure  6).  We  have  therefore  obtained  three  best-fit  linear 
relations  by  the  least-squares  method  to  explain  the  data  from  700  to 
4500  km.  In  the  distance  ranges  of  700  to  1900  km  and  2600  to  4500  km, 
we  have  usee'  all  available  data  for  the  best-fit  relations,  while  between 
2050  and  2550  km  we  have  used  the  earthquake  and  explosion  values  separately, 
he  best-fit  relations  are  given  below 


t = (11. 71+0. 040)+a/(8. 29+0. 002),  700  < A < 1900 

(1) 

Earthquake  Data 

t = (63.37+0.250)+a/(10.44+0.011), 

(2a) 

2050  < A < 2550 

Explosion  Data 

t = (58.43+0.711  )+A/(10. 39+0.033), 

(2b) 

t = (106. 96+0.062)+a/(12. 58+0. 002),  2600  < A < 4500 

(3) 

where  a is  in  kilometers.  The  plots  of  the  above  relations  are  shown  by 

solid  lines  in  Figures  5 and  6.  The  difference  in  apparent  velocities 
in  equations  (2a)  and  (2b)  is  negligible.  This  is  expected  from  the 
parallelism  of  the  two  sets  of  values;  however,  for  reasons  discussed 
previously,  the  relation  in  equation  (2b)  shojld  be  preferred  in  tne 
epicentral  distance  range  in  question. 

Hales  (1972)  summarized  the  results  of  the  seismic  refraction 
studies  in  the  United  States  and  Canada  from  the  Early  Rise  and  NTS 
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explosion  sources.  He  sugcjested  that  betv/aen  600  and  700  I-id,  the  first 
arrival  branch  breaks  off  to  a new  travel-time  branch  with  an  apparent 

velocity  in  the  range  of  8.35  to  8.45  km/sec.  This  branch  extends  up  to 
2100  km  and  yields  an  intercept  time  from  10  to  13  sec.  The  next  travel- 
time branch  extends  from  2100  to  2500  km  with  an  apparent  velocity  of 
approximately  10.0  km 'sec  and  an  intercept  time  of  about  58  sec.  The 
third  branch  (Green  and  Hales,  1968)  beyond  2500  km  shows  an  app?rent 
velocity  of  12.2  km/sec  and  an  intercept  time  of  98.5  sec.  The  breaks 
at  epicentral  distances  2100  and  2500  km  have  been  related  to  the  400 
and  600  km  discontinuities,  respectively.  These  results  indicate  the 
order  of  magnitude  of  the  apparent  P-wave  velocities  and  intercept  times 
which  could  be  expected  on  a regional  basis  in  the  epicentral  distance 
ranges  in  question. 

The  results  presented  from  this  study  permit  a qualitative  comparison 
with  those  cited  above  and  it  can  be  seen  that  the  quantities  in  equations 
(1),  (2b)  and  (3),  which  are  representative  mean  values  over  the  travel 
paths  involved,  show  close  agreement  with  the  regional  values.  Also,  it 
is  important  to  note  that  in  absence  of  the  explosion  data  in  the  epi- 
central distance  range  of  about  2000  and  2600  km,  one  would  be  inclined 
to  relate  the  sharp  changes  in  gradients  at  2100  and  2600  km  with  the 
400  and  600  km  discontinuities,  respectively.  On  the  other  hand,  because 
of  errors  in  location  and  origin  time  for  earthquakes  and  considerable 
uncertainties  in  travel  times  resulting  from  these  sources  for  Aleutian 
earthquakes,  we  would  correlate  the  400  km  discontinuity  with  the  change 
in  travel-time  gradient  at  1900  km  and  the  600  km  discontinuity  with  the 
one  at  2850  km. 
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CONCLUSIONS 


We  cannot  yet  characterize  all  the  crustal  P-wave  phases  from 
local  earthquakes  in  central  Alaska,  due  to  high  scatter  in  data.  We 
hope  to  resolve  this  problem  with  planned  future  explosions  in  the  area 
of  interest.  However,  the  results  obtained  for  the  basal  crustal 
phase  P*  and  P^  could  be  considereo  ob  a representative  mean  for  the 
central  and  southcentral  Alaska. 

The  travel -time  relations  for  the  upper  mantle  P-wave  phases  ob- 
tained from  this  study  are  based  on  source- receiver  configurations  which 
involved  travel  paths  through  different  structural  provinces.  Besides 
this  unavoidable  factor,  inislocations  and  errors  in  origin  time  of 
earthquakes  are  probably  main  sources  of  scatter  in  data.  However,  the 
inclusion  of  data  from  jnderground  nuclear  explosions  for  all  three 
travel-time  branches  provided  a reasonable  basis  to  obtain  the  apparent 
velocities  of  8.29,  10.39  and  12.53  km/sec  for  the  epicentral  distance 
ranges  of  700  to  1900  km,  1900  to  2850  km  and  2850  to  4500  km,  respec- 
tively. They  are  consistent  with  regional  values  as  cited  in  the 
literature.  We  make  no  statement  about  the  presence  or  absence  of  the 
P-LVZ  in  the  upper  mantle  at  this  stage  of  our  studies. 

The  travel  timns  for  the  Aleutian  earthquakes  as  recorded  in  cen- 
tral Alaska  are  systematically  late  by  3 to  4 sec  compared  with  those 
for  underground  nuclear  explosions  on  Amchitka  Island.  By  considering 
krown  results,  we  have  interpreted  this  feature  as  being  due  primarily 
to  mislocations  and  errors  in  origin  time  of  earthquakes  in  that  area; 
the  earthquakes  were  systematically  located  either  deeper  than  true 

focal  depths  or  farther  west  from  true  epicenters  or  a combination  of 
both. 
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The  well  knov/n  400  and  600  km  upper  mantle  discontinuities  have 
been  usually  identified  vvith  the  breaks  in  travel-time  gradients  at 
epicentral  distances  of  2100  and  2600  km,  respectively.  In  tne  absence 
of  available  explosion  data,  one  would  draw  similar  conclusions  as  above 
from  our  results.  However,  weighing  exclusively  in  favor  of  explosion 
data,  we  correlated  the  above  discontinuities  witn  the  breaks  in  travel- 
time gradients  at  1900  and  2850  km.  Nevertheless,  their  exact  positions 
must  be  confirmed  by  incorporating  amplitude  data  besides  travel  times 
and  inverting  the  observed  results  in  terms  of  a realistic  earth  model. 
These  studies  are  in  progress  and  the  results  will  be  reported  in  the 
future. 
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TAIiLi:  I 


REGION 


California 

Off-Coast  of 

British  Columbia 


Southeastern 

Alaska 


Gulf  of  Alaska  and 
Southern  Alaska 


Alaska  Peninsula 


Aleutian 

Islands 


Kurile  Islands 
and 

Kamchatka 


NO.  or 
EVENTS 


12 


5 

39 

10 

5 

22 

3 


rOCAL  DEPTH  (KM) 
9 ^ h jc  14 

12  ^ h £ 33 

18  ^ h jc  25 

21  £ h 126 

20  _<  h ^ 28 

11  ^ h _<  29 
32  _<  h <_  60 

62  _<  h <_  107 

15  _<  h ^ 30 
32  _<  h _<  60 

81  < h < 134 


M/.GNITUDE 


5.1  < m,  < 6.2 

— D — 


4.6  < m.  < 5.9 
— b — 


5.1  ^ ^ 6. 5 


4.0  < m.  < 5.8 
— b — 

majority  > 4.5 
4.8  1 mj^  5.2 


5.0  _<  6.0 

4.6  < m,  < 6.0 
— b — 

majority  nij^  > 5.0 

5.1  m^^  ^ 6.0 

5.1  ^ m^  6.2 
5.0  ^ ^ 6.1 

5.0  ^ < 6.3 
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TABLE  CAPTIONS 


Table  1 


Figure  1 
Figure  2 

Figure  3 


Figure  4. 
Figure  5. 

Figure  6. 


. Ranges  of  magnitude  and  focal  depth  of  earthquakes,  by  region, 
used  in  this  study. 

FIGURE  CAPTIONS 

. Major  tectonic  features  of  central  Alaska. 

. P-wave  travel  time  residuals  for  explosions.  The  short  broken 
lines  indicate  the  depth  and  trend  of  the  Aleutian  deep  seismic 
zone  (Jacob,  1972). 

. a)  Location  of  epicenters  of  earthquakes  used  for  P*,  P^  study. 

b)  P-wave  travel  times  for  earthquakes  for  80  km  < a < 600  km. 

The  locations  and  origin  times  of  earthquakes  were  obtained  using 
the  Jeffreys-Bullen  tables. 

c)  Reduced  P-wave  travel  times  from  earthquakes  for  80  km  < a < 

600  km.  The  locations  and  origin  times  of  earthquakes  were  obtained 
using  the  Jeffreys-Bullen  tables. 

Location  of  epicenters  cf  earthquakes  and  underground  nuclear 
explosions  and  short  period  seismic  stations  used  in  this  study. 
P-wave  travel  times  from  earthquakes  and  explosions  for  700  km 
4500  km.  The  locations  and  origin  times  of  earthquakes 
were  obtained  using  the  Jeffreys-Bullen  tables. 

Reduced  P-wave  travel  times  from  earthquakes  and  explosions  for 
700  km  £A£  4500  km.  The  locations  and  origin  times  of  earth- 
quakes were  obtained  using  the  Jeffreys-Bullen  tables. 
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Major  tectonic  reatures  of  central  Alaska. 


Figure  4.  Location  of  epicenters  of  earthquakes  and  underground  nuclear 

explosions  and  short  period  seismic  stations  used  in  this  study. 


€PtC€NT€R 


Figure  5.  P-wave  travel  times  from  earthquakes  and  explosions  for  700  km 
4500  km.  The  locations  and  origin  times  of  earthquakes 
were  obtained  using  the  Jsffreys-Bullen  tables. 
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2.  SEISMICITY  AND  PLATE  TECTONICS  IN  SOUTH  CENTRAL  ALASKA  I 

By  J.  D.  Vanllormer,  J.  N.  Davies  and  L.  D.  Gedney  I 


ABSTRACT 

Hypocenter  distribution  shows  that  the  Benioff  zone  associated  with 
the  Aleutian  arc  terminates  in  interior  Alaska  some  75  km  north  of  the 
Denali  fault.  There  appears  to  be  a break  in  the  subducting  Pacific 
plate  in  the  Yentna  River-Prince  William  Sound  area  which  separates  tv/o 
seismically  independent  blocks,  similar  to  the  segmented  structure 
reported  for  the  central  Aleutian  arc. 
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INTkODUCTIOtl 


The  role  oC  the  Aleutian  arc  and  the  Denali  fault  in  the  tectonics 
of  the  northeast  I’acific  hasin  has  long  been  a subject  of  study  and 
speculation.  It  is  our  purpose  in  this  report  to  di.scuss  the  term- 
ination of  the  arc  system  in  the  continental  block  and  to  present  some 
new  data  conccrniiuj  the  geometry  of  the  benioff  zone.  First  we  present 
a brief  review  of  past  v/ork,  which  is  intended  only  to  provide  back- 
ground information  and  is  not  meant  to  be  a complete  bibliography. 

The  observations  and  hypotheses  of  some  early  v;orkers  embodied  many 
of  the  elements  of  plate  tectonic  theory,  although  it  was  many  years 
before  the  concept  evolved  into  the  presently  popular  working  hypothesis 
involving  rigid  plates.  Some  twenty  years  ago  Gutenberg  and  Richter 
(1954,  page  32)  commented  that  "The  shallow  shocks  in  the  interior  of 
Alaska  represent  an.  interior  structure,  related  to  the  Pacific  coastal 
arc  as  the  Rochy  Mountains  arc  related  to  t)»e  pacific  coast  ranges. 
Further,  they  observed  that  the  volcanoes  associated  with  the  Aleutian 
Islands  persist  northeastward  onto  the  continent  along  the  Alaska 
Peninsula,  accompanied  by  shocks  of  100  km  depth  and  more.  Benioff  (1954) 
noted  that  earthquake  foci  in  the  Aleuticins  dipped  nortliward  under  the 
arc  and  attributed  the  mechanism  to  a complex  reverse  fault.  These 
fundamental  observations  of  a dipping  seismic  zone  under  the  Aleutian 
arc  with  related  shallow  seismicity  in  the  Alaska  interior  are  still 
valid.  St.  Amand  (1957)  published  a broad  review  and  discussion  of 
Alaska  tectonics  and  applied  the  term  "Denali  Fault  to  the  great 
structural  feature  which  arcs  through  the  Alaska  Range  and  which  had 
been  previously  recognized  in  several  separate  areas  as  a major  fault. 

St.  Amand  also  pointed  out  that  the  Aleutian  trencli  split  just  west  of 
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Kodiak  witlj  the-  decj>watcr  part  branching  toward  the  east  while  the  other 
part  wont  through  Shelikof  Strait,  Cook^lnlet,  and  the  Susitna  depression 
(Figures  1 and  4).  The  implication  of  this  is  that  the  Alaska  Range 
southwest  of  the  "McKinley  Con  " belongs  structurally  to  the  Aleutian 
Arc.  St.  Amand  speculated  that  the  former  trend  of  the  arc  was  through 
Kodiak  Island  and  the  Kenai  Peninsula  to  the  Denali  fault  near  Mt.  Hayes 
while  the  present  trend  is  along  the  Aleutian  Range  to  Mt.  McKinley. 
Counterclockwise  rotation  of  the  northeast  Pacific  basin  about  a pole  in 
the  Gulf  of  Alaska  was  presented  as  a possible  explanation  for  this 
shift  in  trend.  It  is  interesting  to  note  that  relative  movements 
observed  and  implied  by  St.  Amand  in  southcentral  Alaska  are  quite 

similar  to  the  plate  motions  recently  proposed  by  Grow  and  Atwater 
(1970). 

Locating  the  nerthern  terminus  of  the  Benioff  zone  in  Alaska  has 
received  little  attention,  although  several  authors  have  commented  that 
the  intermediate  depth  events  do  not  extend  far  inland.  Tobin  and  Sykes 
(1966)  observed  that  the  intermediate  depth  earthquakes  extended  up  the 
Alaska  Peninsula  250  km  beyond  the  volcanoes  to  the  vicinity  of  Mount 
McKinley,  and  their  map  has  some  events  plotted  just  north  of  the  Denali 
fault.  Tarr  (1970)  also  showed  that  intermediate  depth  events  occurred 
just  north  of  the  Denali  fault,  but  neither  of  the  latter  papers  focused 
much  attention  on  these  events.  Data  presented  in  this  paper  define  the 
northern  terminus  of  the  Benioff  zone. 

Numerous  authors  have  studied  focal  mechanisms  of  Aleutian  arc 
earthquakes  and  the  evidence  for  under-thrusting  of  the  Pacific  Plate  in  a 


•Wc  use  this  phrase  to  denote  the  region  of  the  great 
Alaska  in  the  vicinity  of  Mt.  McKinley. 


structural  bond  in 
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Cross  SL'ctioMS  showing  the  distribution  of  hypoccntors  were  constructed 
by  projecting  nil  the  evfuits  in  n given  rectilinear  volume  onto  a vertical 
plane.  The  volumes  used  were  600  km  in  length,  300  kun  deep  and  25  or  50 
km  thick,  de}»cnding  on  the  density  of  events  in  a given  region.  Earthquakes 
used  in  these  cross  sections  were  limited  to  those  which  had  a computed 
magnitude  of  2.0  or  greater  and  an  rms  travel  time  residual  of  1.0 
second  or  less.  A number  of  events  were  selected  on  tlie  basis  of  magnitude, 
station  distribution  and  dei^th  control  for  processing  by  a modified 
version  of  the  program  HYPO  71  (Lee  and  Lahr,  1972)  and  some  of  the 
results  are  shown  as  section  A7  of  Figure  2.  The  solutions  for  the 
events  which  are  depicted  in  section  A7  are  probably  of  a bettor  quality 
than  many  of  the  ochers.  In  any  case,  the  overall  picture  remains 
unchanged. 

Examination  of  the  cross  sections  reveals  a dipping  zone  of  activity 
chat  is  ofcc;u  quite  well  defined  and  in  some  areas  exhibits  a tendency 
toward  dipping  shallowly  down  to  the  50  km  level  and  then  more  steeply 
below  that  - similar  to  a distinct  change  in  dip  found  in  the  central 
Aleutians  by  Fingdahl  (1971,  1973). 

The  dipping  zone  of  activivy  in  these  cross  sections  is  generally 
80-100  km  thick,  in  rough  agreement  with  the  model  proposed  by  Jacob 
(1972)  and  with  the  observations  of  Engdahl  (1971)  for  the  central 
Aleutian  arc.  Refinement  of  locations  by  modeling  and  ray  tracing  will 
probably  narrow  this  zone  to  some  degree.  Lahr  and  Page  (1972)  and 
Engdahl  (1973)  have  found  the  Benioff  zone  to  be  approximately  10  km 
thick  in  lower  Cook  Inlet  and  the  Central  Aleutians,  respectively. 

Several  sets  of  cross  sections  parallel  to  and  perpendicular  to  the 
Aleutian  Range  and  the  Alaska  Range  (Figures  2A  and  2B)  stiow  that  the 
Benioff  zone  continues  up  the  Aleutian  Range  to  a point  75  km  north  of 
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nort  hv.’er.tcrly  direction  ir;  overwhelming.  As  more  data  have  boon  collected, 
the  simplistic  picture  of  a smootli,  circular  arc  with  a diameter  of  some 
7G0  miles  (Murray,  lO'lS)  has  given  way  to  that  of  a segmental  arc  composed 
of  relatively  short  straight  lines  (V.'ost,  1951)  meeting  at  slight  angles 
to  one  anoUier.  T);ese  blocks  apj)Oar  to  bo  tectonically  independent  and 
they  are  separated  by  ajgiarently  permanent  boundaries  (Jordan  et  al  . , 

19G5,  and  Stauder,  1972).  These  boundaries  have  been  sliown  to  coincide 
with  the  intersections  of  major  transverse  tectonic  features  (Sykes, 

1971)  . There  is  also  a tendency  for  the  main  shock  of  an  eartliquake 
sequence  in  a block  to  occur  near  one  end,  and  for  tlie  aftershock  zone 
to  include  tlie  entire  segment  (Kelleher,  1970;  Sykes,  1971).  St.  Aniand 
(1957)  noted  that  the  maximum  energy  release  in  the  Aleutians  seemed  to 
be  at  junctures  of  the  various  blocks  comj^osing  the  arc. 

DATA  ANALYSIS  AND  INTI- IJ^IXTAT ION 
The  foci  used  in  this  study  were  determined  from  data  collected 
from  stations  of  tlie  seismograph  network  operated  by  the  Geophysical 
Institute,  University  of  Alaska,  and  from  selected  stations  of  the  NOAA 
Palmer  Seismological  Observatory  whicli  wore  recorded  directly  at  the 
University  of  Alaska.  Additional  P-wave  arrival  times  obtained  from 
NOAA/Palmor  for  other  NOAA  stations  not  recorded  at  the  University  of 
Alaska  wore  included  for  events  of  magnitude  larger  than  about  3. 
Hypoconters  (with  the  exception  of  those  shown  in  Fig. 2,  Section  A7,  see 
below)  wore  calculated  by  a least-square  procedure  based  on  P-wave 
tables  compiled  by  Herrin.  (19G8)  . Depths  wore  constrained  to  levels  in 
the  Herrin  tables  and  no  attempt  was  made  to  incorporate  corrections  for 
anomalies  which  may  be  present  duo  to  the  complex  and  largely  unknown 
geologic  structure  of  central  Alaska.  Accuracy  of  epicenters  is  difficult 
to  determine  precisely,  but  all  indications  arc  that  they  are  within  10 
km  of  tJioir  true  positions  (VaiiWormor  ct  al.  , 1973). 


the  Denali  Ciult  where  tli(>  zone  abruptly  teminates  near  the  nortliern 
footliill.  of  the  Alaska  Range  at  about  C4.2h°N  (Figure  2A'j)  . This 
transection  of  the  Denali  fault  by  the  Benioff  zone  casts  doubt  on  the 
interpretation  of  the  Denali  as  a simple  surface  expression  of  tlie 
transform  boundary  between  tlie  Pacific  plate  and  the  continent  analogous 
to  the  San  Andreas  fault  in  California.  Also,  the  extension  of  the 
Denali  fault  system  to  the  west  of  Mt.  McKinley  is  difficult  to  explain 
in  tliis  context.  On  the  otlicr  hand,  it  is  inconceivable  that  such  a 
large,  recently  active  fault  does  not  play  an  important  tectonic  role  in 
the  region  of  the  McKinley  corner.  The  nature  of  this  role  is  not  clear 
at  this  time. 

The  geometry  of  the  dipping  focal  zone  as  determined  from  some  60 
cross  sections  is  illustrated  by  the  contour  map  of  its  upper  surface 
shown  in  Figure  3.  The  Bezioff  zone  is  quite  well  defined  in  this 
manner  except  in  the  vicinity  of  the  Yentna  i.iver  \;hcrc  a relative 
paucity  of  hypoconters  leads  to  nebulous  interpretation  of  cross  sections 
(Figure  2a6) . The  shape  of  the  zone  in  this  region  is  indeterminate 
from  the  hypocenter  data  alone.  North  of  the  Yentna  River  the  zone  dips 
more  steeply  ('''63®)  and  strikes  N45®E  whereas  south  of  the  river  it  dips 
less  steeply  ('\'56®)  and  strikes  more  northerly  (N13®E)  . 

We  interpret  this  anomalous  area  to  be  a break  in  the  structure 
similar  to  the  break  near  Amchitka  Pass  between  the  Rat  Island  block  and 
the  Delarof  Island  block  (Stauder,  1972).  Although  the  various  other 
observations  listed  below  do  not  individually  suggest  this  relationship 
as  strongly  as  the  hyi^ocenter  information,  when  taken  togetlier  they  lend 
support  to  interiircting  this  kink  in  tlio  zone  as  a block  juncture. 
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i;xiir.iinntion  of  ii  tojiA^cji  aj  liic  map  {or  f.ee  1'j.cjui  i-  'I)  £;how;;  that  the 
Yeiitna  Kiver  is  at  a bend  in  Ute  moiiutains  and  that  the  V’idth  of  the 
massif  .soutli  of  the  river  is  approximately  tv;icc  the  width  of  the  massif 
northc'ast  of  the  river.  Gravity  data  for  this  area  (V.’oollard  et  al., 

1960)  indicate  a d«.  cp-r.catcd  change  in  trend  in  this  area,  as  docs  the 
termination  of  the  volcano  chair,  at  Mt.  fJjJurr  100  km  southwest  of  the 
Yentna  River.  Tlierc  is  also  a noticeable  decline  in  seismic  activity 
(riguro  2AG)  within  the  Benioff  zone  in  the  vicinity  of  the  Yentna 
River. 

Aftershocks  of  the  great  1964  Alaska  cart)iquake.  v/erc  confiiicd 
southwest  of  an  extension  of  the  lino  of  the  Yentna  River  anomaly  (Raae 
1968;  Algermisscn  et  al.  , 1972).  This  line,  when  extended  oven  ♦’urthor 
southeast,  coincides  with  the  continental  shelf  break  in  the  Cape  Yakataga 
a’-oa. 

The  initial  .slmck  in  1964  war.  very  near  the  northeast  terninus  of 
the  aftershock  zone.  The  pattern  of  a mai.i  shock  near  one  end  of  a 
block  is  often  observed  in  the  Aleutian  arc  farther  to  the  west  where 
the  tectonic  setting  is  much  sinii^lcr  (Kclleher  et  al.  , 1973)  . The  3 
November  1943  magnitude  7.3  earthquake  also  occurred  near  the  Yentna 
River  lino  and  may  be  further  evidence  of  a fundamental  structural 
break.  It  is  interesting  to  note  that  in  historic  time  no  large  earthquakes 
have  occurred  betv/een  the  Yentna  River  line  and  the  Alaska  Range  to  the 
north  and  northeast  (Davis  and  Echols,  1962). 

Earthquake  magnitudes  determined  by  the  Geophysical  Institute  arc 
generally  higher  for  most  of  Alaska  than  those  reported  by  the  NOAA 
Palmer  Observatory.  Nine  of  eleven  exceptions  to  this  generalization 
for  a randomly  selected  two  month  period  occurred  in  the  region  bounded 
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on  Uie  no;  ,;i  by  tlic  tcniiiimii  of  the  Bciiioff  zone  jncl  on  the  souLli  by  tlie 
Yontna  Kivor  (VanWormcr  et  al.,  1973).  We  liypotlicLiize  tliat  tlie  relatively 
larycr  niacjni tuclen  NOAA  com])utcd  for  this  area  are  duf  to  very  low 
attenuation  of  body  waves  traveling  upward  along  the  dipping  zone  to  the 
KO/Vi  station  at  Palmer.  A similar  situation  of  very  lov;  attenuation  of 
amplitude  was  rcjjorted  for  waves  traveling  up-dip  in  the  Tonga-Kermadec 
arc  by  Oliver  and  Isacks  (1967). 

It  is  clear  from  tlie  above,  observations  of  topography,  gravity,  volcanoes, 
aftershocks,  magnitudes,  and  particularly  hypocenter  distributions  that 
there  is  some  structural  anomaly  at  depth  in  the  Yentna  River  region. 

CONCLUSIONS 

Our  interpretation  of  the  tectonics  in  southcentral  Alask..  s 
schematically  represented  in  Figure  5.  We  propose  that  there  are  two 
blocks  vdiich  arc  similar  in  character  to  the  blocks  composing  the  western 
Aleutian  arc.  Tl>e  boundary  between  these  blocks  extends  from  the  Yentna 
River  through  Prince  William  Sound,  perhaps  as  far  as  the  continental 
shelf  break  near  Yakutat  Bay.  The  northern,  or  McKinley  block,  extends 
northv;est  beyond  the  Denali  fault  to  the  northern  edge  of  the  Alaska 
Range  massif  and  dips  more  steeply  than  the  adjacent  Kcnai  block.  The 
southwest  boundary  of  the  Kenai  block  is  not  presently  known  but  we 
speculate  that  it  may  be  at  the  bend  in  the  volcano  line  in  Kamishak 
Bay,  although  tlie  1964  aftershocks  indicate  the  boundary  may  bo  near  the 
southwest  end  of  Kodiak  Island.  The  evidence  presented  hero  and  by 
Davies  and  Borg  (1973)  indicates  that  the  subduction  zone  in  the  eastern 
Aleutian  Arc  is  tlie  Sho'likof  Strait  - Cook  Inlet  - Susitna  River  topographic 
low  and  is  not  offshore  in  the  trench  area  east  of  Kodiak  Island.  The 
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contiiicnLiil  rnaLfrial  bcLwi;t‘n  Llie  abyssal  plain  and  the  subduction  zone 
is  hypoUiesizcd  Lo  bo  a "raft"  of  material  ridiny  on  the  Pacific  i^late. 
We  speculate  that  the  great  19G4  earthquake,  a shallow  thrust  mecjianism 
(e.g.  Plafker,  1965),  was  a result  of  continued  under thrusting  of  the 
Pacific  plate  beneath  the  rafted  continental  material  of  the  Kenai 
block. 


33 


r 

ackmowli:dg!j:u:i4TS 

vie  express  our  thanks  to  Dr.  D.  H.  Stone  v.’ho  critically  reviewed 
the  manuscript.  This  study  was  funded  by  the  State  of  Alaska,  by  the 
Advanced  Hescarcli  Projects  Ayency  of  the  Department  of  Denfense  (monitored 
by  the  Air  Force  Office  of  Scientific  Research  under  contract  F44620-73-C-0053) , 
and  the  National  Aeronautics  and  Space  Administration  contract  NAS  5-21833. 


WvFEHENCKS 


Algarmi;:scn,  S.  T.  , W.  A.  Rinehart,  R.  VJ.  Sherburne  and  W.  II.  Dillingor,  Jr. 

(1972)  . Prenhocks  and  aftershocks,  in  The  Groat  Alaska  l■:^lrthquakc  of  1964 ; 
Seismology  and  Geodesy,  Vol.  2,  Nat.  Acad.  Sci.,  Wash.  D.  C.,  313-364. 

Benioff,  II.  (1954).  Orogenesis  and  deep  crustal  structure-additional  evidence 
from  seismology,  Gool . Soc.  Am.,  65,  305-400. 

Davies,  J.  and  E.  Derg  (1973).  Crustal  mori^hology  and  plate  tectonics  in 
south  central  Alaska,  Pall.  Seism.  Soc.  Am.  , 63,  673-677. 

Davis,  T.  Neil,  and  C.  Echols  (1962).  A table  of  Alaskan  earthquakes,  1708-1961, 
Geophysical  Research  Report  No.  8,  University  of  Alaska,  Fairbanks,  Alaska. 

Engdahl,  E.  R.  (1971).  Explosion  effects  and  earthquakes  in  the  Amchitka  Island 
region.  Science,  173,  1232-1235. 

Engdahl,  E.  R,  (1973).  Relocation  of  intermediate  depth  earthquakes  in  the 
cc?otral  Aleutians  by  seismic  ray  tracing.  Nature  Phys.  Sci.  , 245,  23-25. 

Grow,  J.  A.  and  T.  Atwater  (1970)  . Mid-Tertiary  tectonic  transition  in  tlie 
Aleutian  Arc,  Geol.  Soc.  Am.,  81,  3715-3722. 

Gutenberg,  B.  and  C.  F.  Richter  (1954).  Seismicity  of  the  Earth  and  Associated 
Phenomena,  Irinceton  University  Press,  Princeton,  New  Jersey,  p.  32. 

Herrin,  F.  (chairman)  (1968).  1968  Seismological  tables  for  P phases, 

Bull.  Seism.  Soc.  Am.,  58,  1196-1211. 

Jacob,  K.  H.  (1972)  . Global  tectonic  imp]  ications  of  anmalous  seismic  P 

traveltimes  from  the  nuclear  explosion  "Longshot",  J.  Geophys.  Res.  , 77, 
2556-2'.:73. 

Jordan,  J.  N.  , R.  A.  Black  and  J.  F.  Lander  (1965)  . Aftershocks  of  the  4 
February  1965  Rat  Island  Earthquake,  Science,  148,  1323-1325. 

Kelleher,  J.  , L.  Sykes,  J.  Oliver  (1973).  Possible  criteria  for  predicting 
earthquake  locations  and  their  application  to  major  plate  boundaries  of 
the  Pacific  and  the  Caribbean,  J.  Geophys.  Res.  , 78,  2547-2585. 


35 


Kclli'lior,  J.  A.  (lyVO).  Spaco-timc:  yeiianiciiy  of  tho  Alaska-A] out j an  .seiumic 

zone,  J.  Geophyr. . , 75,  5745-575C.  j 

I 

I 

Lalir,  J.  C.  and  R.  A.  Pacjo  (1972).  Hypocontral  locations  in  tho  Cook  Inlet  j 

region  of  Alaska,  Abstract  in  1)05; , 53,  1042.  | 

Leo,  W.  U.  K.  , and  J.  C.  Lalir  (1972).  I!YP071:  a computer  program  for  determining 

hyjiocentcr,  magnitude,  and  first  motion  pattern  of  local  earthquakes, 

I 

US  Geolociical  Ouivcy,  0]:oii-filo  Report.  100  p.  | 

Murray,  H.  W.  (1945)  . Profiles  of  tho  Aleutian  Trench,  Gool.  Soc.  Am.  , 5G,  I 

757-782. 

Oliver,  J.  and  B.  Isacks  (1967).  Deep  Earthquake  zones,  anomalous  structures  in 
the  upper  mantle,  and  tho  lithosphere,  J.  Goophys.  Res. , 72,  4259-4275. 

Pago,  R.  (1968)  . Aftershocks  and  microaf torsliocks  of  tho  groat  Alaska 
earthquake  of  1964,  Bull.  Seism.  Soc.  Am.  , 58,  1131-1168. 

Plafker,  G.  (1965).  Tectonic  deformation  associated  with  tho  1964  Alaska 
earthquake.  Science,  148,  1675-1687. 

St.  Amand,  P.  (1957)  . Geological  and  geophysical  synthesis  of  tho  tectonics  of 

portions  of  British  cclurabia,  the  Yukon  Terriory,  and  Alaska,  C-eol.  Soc.  Am., 

68,  1343-1370. 

Staudor,  W.  (1972).  Fault  motion  and  spatially  bounded  cliaractcr  of  earthquakes 
in  Amchitka  Pass  and  tho  Dclarof  Islands,  J.  Goophys.  Res. , 77.  2072-2080. 

Sykes,  L.  R.  (1971).  /vftorshock  zones  of  great  earthquakes,  seismicity  yaps, 
and  earthquake  prediction  for  Alaska  and  tho  Aleutians,  J.  Goophys.  Res., 

76,  8021-8041. 

Tarr,  A.  C.  (1970).  Now  Maps  of  Polar  Seismicity,  Bull.  Seism.  Soc,  Am.  , 60, 

1745-1757. 

Tobin,  D.  G.  and  L.  R.  Sykes  (1966).  Relations  of  hypocentors  of  earthquakes  to 
the  geology  of  Alaska,  J.  Goopiiys.  Res.  , 71,  1659-1668. 


36 


V..,«t,rm,  r,  D.  , J.  ftwica  „„d  L.  Cc^cl.ioy  (1973).  Central  Alaska  oarthquakas  <luri„.j 

19)2,  ££j.._Kar>l.  . llAC_lj^2<  Ge-oi.liyBical  In.-.titutc,  University  of  Alaska, 
Fairbanks,  Alaska. 

West,  S.  S.  (1951).  Major  shear  fractures  of  Alaska  and  adjacent  regions, 

Trans.  Ain.  Geophys.  Union.  32,  Gl-96. 

Woollard,  G.  P.,  N.  A.  OsLenso,  E.  Thiel,  and  W.  E.  Bonini  (1960).  Gravity 

anomalies,  crustal  structure,  and  geology  in  Alaska,  J.  Geophys.  Res. . G5, 
1021-1037. 


37 


Figure  Captions 

Figure  1.  Location  map  of  Alaska  showing  major  features  referenced  in 
text. 

Figure  2A.  Cross  sections  through  the  Benioff  zone.  Note  that  the 
vertical  exaggeration  is  2:1  in  aJl  sections.  Sections  A1-A6  arc  50  km 
thick  perpendicular  to  the  plane  of  the  figure.  Section  hi,  500  km 
thick  perpendicular  to  the  plane  of  the  figure,  illustrates  the  clarity 
of  the  dipjiing  zone  wiicn  strict  quality  control  is  applied  to  the  data. 

Locations  and  Icok  directions  of  tlie  crocs  sections  are  shown  in  the 
accomi)anying  maps.  Note  the  abrupt  termination  of  the  Benioff  zone  to 
the  left  (north)  in  section  A6. 

figure  2B.  Cross  sections  of  the  nortlicrn  part  of  the  Benioff  zone. 

These  sections  are  25  kin  thick  perpendicular  to  the  plane  of  the  figure 
and  are  nearly  along  the  strike  of  the  zone.  Locations  and  look  directions 
are  shown  in  the  accompanying  map.  Sections  A1-A6  and  B1-B8  are  representative 
of  the  more  than  60  cross  sections  analyzed  in  this  study. 

Figure  3.  Contours  of  the  upper  surface  of  the  Benioff  zone  as  detoi-mined 
in  this  study.  Contour  interval  is  50  km.  Southwest  termination  of  the  contours 
represents  only  the  arbitrary  limit  of  hypocenter  selection.  Northeast 
termination  of  contours  represents  real  termination  of  the  Benioff  zone.  The  star  is 
the  epicenter  of  the  1964  earthquake  and  the  solid  circle  icthc  epicenter 
of  the  1943  m^gritude  7.3  c.arthquake  (as  rclcca*-cd  by  Sykes,  1971). 

Figure  4.  Mosaic  of  19  Earth  Resources  Technology  Satcilite  (ERTS-1)  images. 

Images  were  obtained  in  early  November,  1972.  The  sun  is  illuminating 
the  area  from  the  lower  right.  The  Denali  fault  arcs  across  the  scene  from 
upper  right  to  center  left. 

Figure  5.  Diagrammatic  representation  of  the  subducting  Pacific  plate  in 
south  central  Alaska  as  interpreted  in  tliis  study.  Circle  and  star  are  as  in 
Figure  3.  As  stated  in  the  text,  the  precise  nature  of  the  Benioff  zone  at 
the  Ycntna  lineament  is  not  clearly  defined.  Thus  the  depth  contours  in 
the  figure  arc  broken  at  the  bend  in  the  zone. 
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Figure  2A.  Cross  sections  through  the  Benioff  7.cnc.  Note  that  the 
vertical  exaggeration  is  2:1  in  all  sections.  Sections  A1-A6  arc  50  km 
thick  perpendicular  to  the  plcine  of  the  figure.  Section  A7,  500  km 
thick  iierjiendicular  to  the  plane  of  the  figure,  illustrates  the  clarity 
of  the  dii'ping  none  when  strict  quality  control  is  applied  to  the  data. 
Locations  and  look  directions  of  the  cross  section.s  are  shown  in  the 
accoiijAinying  maps.  Note  the  abrupt  termination  of  the  Benioff  zone  to 
the  left  (nortli)  in  section  A6. 
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Figure  3.  Contours  of  the  upper  surface  of  the  Benioff  zone  as  determined 
in  this  study.  Contour  interval  is  50  km.  Soutliwcst  termination  of  the  contours 
represents  only  the  arbitrary  limit  of  iiypocenter  selection.  Northeast 
termination  of  contours  rc  prei-.cnts  real  termination  of  the  Benioff  zone.  The  star  is 
the  epicenter  of  the  lO&'l  earthquake  and  the  solid  circle  isthe  epicenter 
of  the  1043  magnitude  7.3  earthquake  (as  relocated  by  Sykes,  1971). 


Figure  A.  Mor.aic  of  19  FarLli  Rcsourccrj  Technology  Satellite  (FRT 
Images  v/cre  obtained  in  early  Novemlx^r,  1972.  The  sun  is  illumin 
the  area  from  the  lower  right.  The  Denali  fault  arc:  across  the 
upper  rJg'at  to  center  lett. 


Fiqurc  5.  Diagrammatic  representation  of  the  subducting  Pacific  plate  in 
south  central  Alaska  as  interpreted  in  this  study.  Circle  and  star 
Figure  3.  As  slated  in  the  text,  the  precise  nature  of  the  Bonioff  zone  at 
the  Yentna  lineament  is  not  clearly  defined.  Thus  the  depth  contours  in 
the  figure  arc  broken  at  the  bend  in  the  zone. 


3.  Source  Mechanics  for  Small  Magnitude  Earthguakes  in  Central  Alaska 

by  B.  Bliattacharya  and  N.  N.  Biswas 

ABSTRACT 

P-wave  first-motion  for  41  selected  earthquakes  (3.0^m^^5.5)  located  in 
central  and  southcentral  Alaska  have  been  studied.  The  events  with  focal  depths 
less  than  60  km  show  either  strike-slip  or  normal  faulting  focal  mechanism. 

The  events  with  focal  depths  greather  than  the  above  value  show  nearly  pure 
dip-slip  motions.  The  solutions  in  general  show  a northwest-southeast 
orientation  of  the  pressure  axes,  a feature  in  good  agreement  with  that 
predicted  from  North  Pacific  plate  tectonics. 
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INTRODUCTION 


The  Anchorage  earthquake  of  1964  led  to  more  research  aimed  at 
understanding  the  seisinotectonic  setting  of  the  Alaska-Aleutian  region. 

Some  of  the  important  results  revealed  by  the  ensuing  studies  and  which 
are  pertinent  to  the  present  investigation  have  been  discussed  by  Stauder 
(1968a»b).  From  the  first-motion  solutions,  Stauder  demonstrated  a 
close  correlation  between  the  nature  of  motions  on  active  faults  along 
the  coastal  belt  of  central  Alaska  and  Aleutian  regions  and  those  predicted 
from  the  hypothesis  of  global  tectonics  for  the  above  regions. 

Relationships  between  fault-plane  solutions  and  regional  tectonics 
in  central  Alaska  have  been  investigated  by  Gedney  and  Berg  (1969a, b). 
Gedney  (1970)  integrated  earlier  fault-plane  solutions  with  those  from 
small  magnitude  earthquakes  located  from  Kodiak  Island  to  the  Alaska 
Range  along  a NNE-SSW  line  through  Cook  Inlet.  For  a number  of  events, 
he  obtained  a northwest-southeast  ttending  pressure  axis,  a feature 
which  is  in  agreement  with  the  results  obtained  for  the  coastal  belt  of 
central  Alaska  by  Stauder  and  Bollinger  (1966)  and  Stauder  (1968b). 

On  the  north  of  the  Alaska  Range,  fault  plane  solutions  for  two 
shallow  moderate  size  earthquakes  and  their  aftershock  sequences  have 
been  discussed  by  Jordan  et  al.  (1968),  Gedney  and  Berg  (1969)  and 
Gedney  (1970).  The  first  event  was  located  about  100  miles  west  of 
Fairbanks  and  showed  a left-lateral  movement  along  a north-south  trending 
steeply-dipping  fault.  The  second  event  was  located  just  south  of 
Fairbanks  and  occurred  along  a northwest-southeast  trending  fault;  the 

motion  had  equal  components  of  right-lateral  strike-slip  and  normal 
faulting. 


Bisv/as  (1973),  from  P-wavo  travel  time  residual  studies,  inferred 
the  inland  continuation  of  the  Aleutian-Alaska  Peninsula  underthrust 
lithosplieric  plate  under  southcentral  and  central  Alaska.  As  incor- 
porated in  Section  2 of  this  report,  Van  Wonner  et  al.  (1974)  have 
demonstrated  the  termination  of  the  above  feature  some  75  km  north  of 
the  Denali  fault.  Lahr  (1974),  in  addition  to  obtaining  similar  re- 
sults, has  related  the  fault  plane  solutions  in  Cook  Inlet  to  the 
geometry  of  the  plate  from  improved  seismic  coverage  around  the  above 
area.  He  obtained  a N28°E  striking  Benioff  zone  with  a dip  of  about 
45°-50°;  the  plate  was  shown  to  reach  a depth  of  about  150  km  from  the 
surface.  Lahr  showed  also  that  the  earthquakes  at  intermediate  depths 
(>  60km)  are  characterized  by  nearly  pure  dip-slip  motions  with  the 
tension  axis  oriented  downdip  along  the  seismic  zone. 

In  the  present  study,  we  have  examined  the  source  mechanisms  of 
small  magnitude  earthquakes  in  central  and  southcentral  Alaska  with 
respect  to  the  current  knowledge  of  the  North  Pacific  plate  tectonics. 

DATA 

We  have  analyzed  short  period  records  for  the  entire  year  of  1973 
obtained  by  the  central  Alaskan  seismic  networks  (Geophysical  Institute 
and  NOAA).  The  correct  polarity  of  the  seismograph  at  each  station  was 
checked  by  studying  the  sense  of  motion  from  a presumed  underground  explosion 
(12  September  1973,  = 6.8)  in  Novaya  Zamlaya. 

Out  of  the  total  of  108  events  located  in  the  area  of  interest,  P- 
wave  first  motions  for  41  events  could  be  read  at  most  of  the  stations 
with  certainty.  The  locations  and  origin  times  for  these  events  were 
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noted  from  NOAA's  Monthly  Bulletin  (PDE)  and  are  listed  in  Tables  1 and 
2;  the  magnitudes  of  the  earthquakes  vary  from  3.0  to  5.5.  The  locations 
of  the  stations  and  epicenters  of  the  events  considered  are  shown  in 
Figure  1. 

To  determine  the  emergence  angle  at  the  focus  for  each  earthquake, 
a plane  layered  crust-upper  mantle  P-wave  velocity  structure  (Table  4) 
has  been  constructed  from  the  results  discussed  in  Section  1 of  this 
report.  Since  we  cannot  characterize  the  P^  crustal  phase  yet,  P- 
wave  velocity  for  the  first  crustal  layer  has  been  taken  from  Hales  and 
Asada  (1966).  Although  this  value  for  the  P^  phase  is  appropriate  for 
southcentral  and  southeastern  Alaska,  we  do  not  anticipate  any  significant 
errors  in  using  the  above  P^  value  for  all  of  central  Alaska  in  computing 
the  angle  of  emergence. 

The  computations  of  the  emergence  angles  have  been  performed  by  the 
HYPO-71  (Lee  and  Lahr,  1972)  computer  program.  For  most  of  the  events, 
the  combination  of  short  epicentral  distance  and  intermediate  focal 
depth  yielded  the  first  arrivals  as  direct  waves.  Due  to  this  fact,  we 
have  represented  the  first-motion  dilatations  and  compressions  on  an 
equal  area  projection  (10  cm  Schmidt  net)  of  the  upper  hemisphere.  For 
those  cases  where  the  first  arrivals  were  refracted,  the  representation 
on  the  upper  hemisphere  has  been  retained  for  consistency. 

RESULTS 

The  representation  on  the  stereo-net  of  the  distribution  of 
compressions  and  dilatations  of  first-motions  for  each  event  was  visually 
best-fitted  to  a nodal  plane.  A second  nodal  plane,  orthogonal  to  the 
first,  was  then  drawn  as  a best-fit  through  the  data  points.  This 
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procedure  was  applied  to  the  data  for  all  the  41  earthquake's;  the 
results  have  been  divided  into  two  groups  on  the  basis  of  focal  depth 
and  the  nature  of  the  solution,  that  is,  dip-slip  or  other  type,  the 
discussion  of  which  follows. 

The  earthquakes  listed  in  Table  1 have  hypocenters  located  either 
in  the  crust  or  10  to  15  km  deeper  than  the  average  depth  of  the  M- 
discontinuity  45  km)  in  central  Alaska  (Table  4).  These  have  been 
placed  in  one  group.  Their  first  motion  solutions  are  shown  in  Figures 
2 and  3,  which  show  either  strike-slip  or  normal  faultings.  The  hypo- 
centers  of  the  earthquakes  listed  in  Table  2 vary  in  depth  from  64  to  171 
km,  with  one  exception  (event  No.  20),  for  which  it  is  32  km. 

However,  for  the  majority  of  the  events,  the  focal  depths  exceed  100  km. 
The  first-motion  solutions  of  all  these  events  show  nearly  pure  dip-slip 
motions  as  shown  in  Figures  4 and  5.  These  solutions  have  been  placed 
in  the  second  group. 

In  Tables  1 and  2 we  have  also  shown  the  values  of  the  dip  direction, 
dip  and  slip  angle  (angle  measured  in  the  nodal  plane  between  the  hori- 
zontal and  the  di»-ection  of  motion,  if  the  nodal  plane  in  question  be 
considered  as  the  fault-plane)  for  each  of  the  nodal  planes  and  the 
trend  and  the  plunge  of  the  B,  P and  T axes.  The  B axis  is  the  null 
axis,  i.e.  the  intersection  of  the  two  nodal  planes.  The  P and  T axes 
are  the  axes  of  greatest  and  least  compressive  stresses,  respectively. 

Of  the  two  nodal  planes,  plane  1 refers  to  the  one  which  has  relatively 
greater  dip  value. 

In  Figures  6 and  7 the  first-motion  solutions  of  the  first  and 
second  groups  of  events,  respectively,  have  been  shown  in  relation  to 
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the  major  structural  trends  in  central  Alaska.  These  include  seven 
prominent  arcuate  fault  systems  with  right  lateral  motions,  except  the 
faults  shown  by  barbed  lines  along  the  coast  which  are  thrust  type. 

The  shaded  and  non-shaded  sections  in  the  mechanism  diagrams  represent 
the  dilatational  and  compress ■'onal  areas,  respectively.  The  trend  of 
the  pressure  axis  for  each  case  is  shown  by  a solid  arrow. 

Six  of  the  events  (Nos.  1,8,  10,  11,  14  and  18)  of  the  first  group 
represent  strike-slip  motion  (Figure  6);  one  of  the  nodal  plane  strikes 
approximately  in  the  east-west  direction,  while  the  other  in  the  north- 
south  direction  and  both  of  them  dip  steeply  by  more  than  80°.  The 
pressure  axes  are  nearly  horizontal  (1°  < plunge  < 8°,  Table  1)  and  are 
oriented  either  in  the  southeast-northwest  (Nos.  8,  10,  11,  18)  or  in  the 
northeast-southwest  (Nos.  1,  14)  direction. 

Three  (Nos.  1,  11  and  18)  of  the  strike-slip  events  are  located  in 
the  vicinity  of  the  faults  marked  H,  C and  G (Figure  6).  On  the  basis 
of  the  orientation  of  the  fault  traces,  v/e  interpreted  the  east-west 
striking  nodal  plane  as  the  fault-plane  for  these  events.  The  remaining 
three  events  (Nos.  8,  10  and  14)  are  located  in  valleys  having  thick 
sedimentary  deposits  with  no  exposed  fault  traces.  Thus  we  have  not 
identified  which  of  the  two  nodal  plane^  is  the  fault  plane  for  these 
cases. 

The  first-motion  solutions  for  the  remaining  twelve  earthquakes  of 
the  first  group  represent  normal  faulting.  Except  for  Nos.  2 and  7,  the 
events  show  the  dip  of  the  two  nodal  planes  in  the  range  of  30°  to  60°. 
For  the  two  former  events,  the  dip  increases  to  69°  to  71°.  The 
plunge  of  the  tension  axis  has  been  found  to  vary  from  0°  to  28°. 
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Moreover,  for  six  of  the  events  (Nos.  4,  9,  12,  13,  15,  17)  the  two  nodal 
planes  strike  approximately  in  the  northwest-southeast  direction  whereas 
four  of  them  (Nos.  2,  3,  5 and  6)  show  about  northeast-southwest  and  two 
(Nos.  7 and  16)  show  approximately  east-west  striking  nodal  planes. 

The  earthquakes  of  the  second  group  are  all  located  in  southcentral 
Alaska  (Figure  7).  The  first-motion  solutions  for  these  events  correspond 
to  underthrust  type  of  motion,  that  is,  the  foot  wall  is  sliding  downward 
with  respect  to  the  hanging  wall.  Also,  the  direction  of  the  dip  of  one 
of  the  nodal  planes  varies  from  west  to  northwest  while  that  of  the  other 
varies  from  east  to  southeast.  Exception  to  these  general  dip  directions 
were  obtained  for  Nos.  5,  6,  13,  17,  21  and  23.  For  these  cases,  one  of 
the  nodal  planes  dips  either  in  the  south  or  southwest  direction  while 
the  other  in  the  north  and  northeast  direction.  The  trend  of  the  nodal 
planes  for  event  No.  8 does  not  belong  to  either  of  the  above  two  types; 
one  of  the  nodal  planes  dips  in  the  NNE  direction,  the  other  in  the  SSW 
direction. 

The  values  of  the  dip  for  both  the  nodal  planes  range  from  30°  to 
60°  except  for  events  No.  1,  16,  22,  for  which  it  is  less  than  30°. 

Also  the  plunge  of  the  tension  axis  varies  from  68°  to  88°,  except  for 
events  No.  1 and  8.  For  these  tv/o  cases,  the  values  are  53°  and 
46°,  respectively.  The  plunge  of  the  pressure  axis  varies  from  b°  to 
22°  except  for  event  No.  1,  for  which  it  is  37°. 

DISCUSSION 

Although  the  extension  of  the  North  Pacific  lithospheric  plate 
across  the  Denali  fault  to  central  Alaska  up  to  at  least  63°N  has  been 
suggested  (Van  Wormer  et  al.,  1974;  Lahr,  1974),  the  relationships  between 
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the  sources  of  high  seismicity,  particularly  from  the  Denali  fault  to  about 
66°N,  and  the  plate  motion  are  difficult  to  resolve  due  to  the  complicated 
structural  setting  of  the  above  section  of  central  Alaska  (Figure  6). 
However,  to  derive  a regional  pattern  from  tfie  individual  first-motion 
solutions,  we  have  prepared  composite  diagrams  for  the  strike-slip  and 
normal  faults  separately  for  events  located  north  of  the  Denali  fault. 
Figure  8a  shows  such  a solution  for  the  strike-slip  faults  revealed  by 
event  Nos.  8,  10,  11,  14  and  18. 

In  the  above  figure,  two  nodal  planes,  one  striking  approximately 
east-west  and  the  other  north-south,  give  a visual  best-fit  to  the  data. 

If  the  former  nodal  plane  is  chosen  as  the  regional  orientation  of  the 
strike-slip  faults,  then  the  solution  v;ill  imply  a right-lateral  motion 
on  a nearly  vertical  fault  with  pressure  axis  oriented  in  the  northwest- 
southeast  direction. 

Figure  8b  shows  the  composite  solution  of  normal  faults  for  earth- 
quakes Nos.  4,  9,  12,  13  and  15.  Data  for  tv;o  of  the  events  (Nos.  5 and 
7)  have  been  excluded  from  the  above  figure  due  to  grossly  different 
strike  directions  of  their  nodal  planes  compared  to  others.  This  does 
not  imply  that  the  first  motion  solutions  for  these  two  events  are 
inaccurate,  rather  the  orientation  of  their  fault  planes  does  not  follow 
the  regional  trend,  a feature  not  uncommon  for  a complicated  active 
tectonic  area. 

The  above  composite  solution  (Figure  8b)  shows  a steeply  dipping 
northwest-southeast  orientation  of  the  pressure  axis,  which  is  in  close 
agreement  with  that  obtained  for  the  strike-slip  faults.  The  visual 
best-fit  of  the  data  with  two  nodal  planes  has  uncertainties  which  are 
equivalent  to  errors  from  1 5°  to  + 10°  in  the  dips  and  strikes  of  the 
nodal  planes. 
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The  first-motion  solutions  of  the  intermediate  depth  earthquakes  in 
southcentral  Alaska  have  a general  characteristic  of  underthrusting, 
i.e.  nearly  pure  dip-slip  motions,  as  mentioned  before.  In  addition, 
if  the  northwest  dipping  nodal  plane  is  interpreted  as  the  fault  plane, 
the  focal  model  for  most  of  the  events  will  correspond  to  an  average 
plunge  of  75°  for  the  tension  axis.  These  features  imply  a downward 
sliding  of  the  foot  v/all  with  respect  to  the  hanging  wall,  a result  in 
close  agreement  with  those  obtained  by  Lahr  (1974)  in  the  Cook  Inlet 
area. 

CONCLUSIONS 

In  view  of  the  current  knowledge  of  the  North  Pacific  plate  tec- 
tonics in  central  Alaska,  we  may  infer  that  the  shallow  earthquakes 
occurring  in  the  above  area  due  to  a direct  consequence  of  the  bending 
of  the  underthrusting  lithospheric  plate,  most  likely  will  be  associated 
with  normal  faulting.  If  the  stored  strain  energy  is  released  on  the 
pre-existing  fauUs,  the  above  pattern  may  change  to  strike-slip  motions. 

On  the  other  hand,  the  source  mechanism  of  earthquakes  with  intermediate 

focal  depth  (>  6t)km),  resulting  from  the  sinking  of  the  lithospheric 

plate  under  its  own  weight  (Le  Pichon  et  al.,  1973)  will  be  associated  with  nearly 

pure  underthrusting  motions.  This  implies  dov<ndip  orientation  of 

the  tension  axes  along  the  seismic  zone.  If  the  rupture  plane  parallels 

the  dip  of  the  seismic  zone,  the  tension  axes  will  be  at  an  angle  ('\-  45°) 

to  the  plate  boundary.  For  small  magnitude  earthquakes,  whether 

the  rupture  plane  will  always  be  parallel  to  the  plate  boundary  or  not 

seems  to  be  an  open  question. 

The  results  of  the  present  study  for  earthquakes  with  focal  depth 
less  than  60  km  in  central  Alaska  are  in  general  agreement  with  the 
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observations  for  tlie  causes  of  shallow  events  made  above.  However,  we 
obtained  a number  of  solutions  in  the  above  group  which  do  not  agree 
with  the  major  tectonic  pattern  of  the  area.  Further  studies  of  source 
mechanism  for  events  located  in  and  around  those  which  yielded  incon- 
sistent solutions  are  needed  for  a complete  understanding  of  the  nature 
of  the  stored  strain  energy  release  in  the  area  of  interest.  Such 
studies  are  particularly  important  for  areas  which  are  seismically 
active  and  at  the  same  time,  lie  outside  the  northern  limit  of  the 

underthrust  lithospheric  plate  in  interior  Alaska.  • 

The  earthquakes  with  focal  depths  greater  than  60  km  are  all  located  I 

in  southcentral  Alaska  where  the  North  Pacific  plate  tectonics  seem  to  \ 

i 

control  the  seismic  activities.  The  focal  mechanism  of  the  above  events 
show  good  correlation  with  the  features  discussed  above  for  an 
unde*'thru';<’  ;il^te  sinking  by  its  own  weight. 
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Locations  (Lat.t  Lonj.),  origin  tines  and  deti'ls  of  the  P-wave  first-notion 
solutions  for  23  earthquakes  (focal  depth  > 6J  km)  considered  In  this  study. 
Event  Date  Origin  Time  Plane  1 Plane  2 B 

•<unber  1973  Hr  Kin  Sec 
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• 'ons  (Lat..  Long.)  of  the  ^hort  period  seismic  stations  in 

Alaska  used  for  this  study. 


central 


Serial 

No. 


Station 

Name 


1. 

Arctic  Valley 

PMS 

2. 

Black  Rapids 

BLR 

3. 

Clear  Creek  Butte 

CCB 

4. 

Engineer  Hill 

ENG 

5. 

Fort  Yukon 

FYU 

6. 

Gilmore 

GLM 

7. 

Gilmore  Creek 

GIL 

8. 

Granite  Mountain 

GMA 

9. 

Hepp 

HPP 

10. 

Houston 

PWA 

11. 

Indian  Mountain 

iHA 

12. 

Kodiak 

KDC 

13. 

Levy 

LVY 

14. 

McKinley 

MCK 

15. 

Mercer 

MCR 

16. 

Nenana 

NEA 

17. 

Palmer 

PMR 

18. 

Paxon 

PAX 

19. 

Remote 

RON 

20. 

Sheep  Mountain 

SCM 

21. 

Sparrevohn 

SVW 

22. 

Tanana 

TNN 

23. 

Tatal ina 

HA 

24. 

Tol sona 

TOA 

Lat. 

(°N) 

Long. 

(°w) 

Elevation 

(Meter) 

61.2 

149.5 

716 

63.5 

145.8 

809 

64.6 

147.8 

183 

64.7 

147.0 

220 

66.5 

145.2 

137 

64.9 

147.4 

722 

64.9 

147.5 

350 

65.4 

161.2 

858 

64.8 

147.9 

170 

61 .6 

149.8 

137 

66.0 

153.6 

1233 

57.7 

152.5 

13 

64.2 

63.7 

63.9 

64.5 

61.6 

62.9 

62.7 

61.8 

61.1 

65.2 

62.9 

62.1 


149.2 

148.9 

149.0 

149.0 

149.1 

145.5 

150.2 

147.3 

155.6 

151.9 

156.0 

146.2 
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TABLE  4 

Layered  crust-upper  mantle  model  used  for  this  study. 


Depth  to  the  Interface 
(km) 

P-wave  Velocity 
{ km/sec) 

0.0 

5.84 

17.6 

7.05 

4b. U 

8.07 

75.7 

8.29 

300.7 

10.39 

544.8 

12.58 
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TABLE  CAPTIO.JS 


Table  1.  Locations  (Lat. , Long.)*  origin  times  and  details  of  the  P-wave 
first-motion  solutions  for  18  earthquakes  (focal  depth  < 60  km) 
considered  in  this  study. 

Table  2.  Locations  (Lat.,  Long.),  origin  times  and  details  of  the  P-wave 
first-motion  solutions  for  23  earthquakes  (focal  depth  > 60  km) 
considered  in  this  study. 

Table  3.  Locations  (Lat.,  Long.)  of  the  short  period  seismic  stations  in 
central  Alaska  used  for  this  study. 

Table  4.  Layered  crust-upper  mantle  model  used  for  this  study. 


riGURE  CAPTIOUS 


Figure  1 . 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


Figure  8a. 


Figure  8b. 


Locations  of  seismic  stations  and  epicenters  of  earthquakes 
considered  for  the  present  study. 

P-v/ave  first-motion  solutions  for  earthquakes  (No.  1-9)  with 

focal  depth  less  than  60  km.  Crosses  (x)  indicate  uncertain 
data  points. 

P-wave  first-motion  solutions  for  earthquakes  (No.  10-18)  with 
focal  depth  less  than  60  km. 

P-wave  first-motion  solutions  for  earthquakes  (No.  1-12)  with 
focal  depths  greater  than  60  km. 

P-wave  first-motion  solutions  for  earthquakes  (No.  13-23)  with 
focal  depth  greater  than  60  km. 

Relationships  of  the  major  structural  trends  in  Alaska  and 
first-motion  solutions  of  the  18  earthquakes  listed  in  Table 

I.  Shaded  areas  indicate  dilatation  and  non-shaded  areas 
compression.  Solid  arrov/s  indicate  the  orientation  of  the 
pressure  axis. 

Relationships  of  the  major  structural  trends  in  Alaska  and 
first-motion  solutions  for  the  23  earthquakes  listed  in 
Table  2. 

Composite  first-motion  solution  for  earthquakes  (No.  8,  10, 

II,  14  and  18) . 

Composite  first-motion  solution  for  earthquakes  (No.  4,  9, 

12,  13  and  15). 
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A SEISMIC  STATION 
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Locations  of  seismic  stations  and  epicenters  of  earthquakes 
considered  for  the  present  study. 
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Figure  5.  P-wave  first-motion  solutions  for  earthquakes  (No.  13-23 
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Relationships  of  the  major  structural  trends  in  Alaska  and 
first-motion  solutions  of  the  18  earthquakes  listed  in  Tabl 
1.  Shaded  areas  indicate  dilatation  and  non-shaded  areas 
compression.  Solid  arrows  indicate  the  orientation  of  the 
pressure  axis. 


Relationships  of  the  major  structural  trends  in  Alaska  and 
first-motion  solutions  for  the  23  earthquakes  listed  in 
Table  2. 
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4.  STUDY  or  r:s-M^  RElATIOf^S  FROM  CENTRAL  ALASKAN  DATA 

by 

N.  N.  Bisv/as  and  B.  Bhattacharya 
ABSTRACT 

Short  and  long  period  seismograms  obtained  by  the  standard  College 
Observatory  in  central  Alaska  from  1968  through  1972  have  been  analyzed  from 
u-v-erground  explosions  and  earthquakes  located  in  and  around  the  areas 
of  explosions  for  two  regions:  Cal ifornia-Nevada  and  the  Aleutian  Island 
chain.  In  the  Ms-mj^  relations,  the  results  for  explosions  showed  clear 
separation  from  those  for  earthquakes. 
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INTRODUCTION 

Many  authors  (Licbermann  and  Pomeroy,  1969;  Evernden  et  a1.,  1971; 
Ward  and  Toksoz,  1971)  have  shown  from  theoretical  and  experimental 
results  that  a significant  difference  in  the  surface  wave  (M^)  versus  body 
wave  (m^)  magnitude  exists  between  earthquakes  and  underground  explosions. 
Also,  regional  variations  of  M^-m^  relationships  for  earthquakes  have 

been  reported  from  time  to  time.  The  results  have  been  summarized  by 
Aki  (1972). 

Representing  an  earthauake  source  as  a displacement  discontinuity 
across  a plane,  Aki  (1967)  showed  that  for  an  infinite,  isotropic  and 
homogeneous  medium,  the  far-field  displacement  spectra  can  be  expressed 
as  a convolution  of  two  terms.  The  first  term  constitutes  the  wave  propa- 
gation  properties  and  does  not  include  the  fault  motion  parameters  except 
the  direction  of  fault  motion.  The  second  term  includes  the  spatial 
and  temporal  dependence  of  f.e  source.  Also,  Aki  showed  that  the 
source  term  decreased  inversely  as  beyond  a corner  frequency,  where 
w is  the  angular  frequency.  Aki  identified  this  as  the  w^-model. 

By  comparing  such  parameters  as  fault  length,  stress  drop,  seismic 
moment,  and  spectral  ratio  for  periods  greater  than  10  sec  and 
relations  from  observed  data  with  those  derived  theoretically  according 

to  the  oj^-model,  Aki  (1972)  showed  close  agreement  between  the  two  sets 
of  results  for  earthquakes  of  M^>6. 

For  M^<6,  results  of  several  studies  indicated  disagreement  with 
those  from  u^^-model  particularly  for  the  M^-m^^  relation.  Moreover,  for 
two  nearby  earthquakes  in  the  Aleutian  area  of  similar  source  mechanism 
but  of  magnitude  7 and  6,  the  corner  frequencies  were  found  in  disagreement 
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vnth  the  j -iiiodpl.  These  results  have  been  reviewed  by  Aki  (1972). 

The  present  study  v;as  initially  proposed  to  investigate:  (1)  the 

sources  of  the  discrepancy  between  the  observed  and  predicted  results 
for  the  u -model,  particularly  for  the  Aleutian  earthquakes  of  M^<6  and 
(2)  regionalized  i"olations  for  earthquakes  and  underground  nuclear 

explosions.  In  these  investigations,  v/e  emphasized  the  use  of  the  high  quality 
local  ALPA  data  for  the  surface  wave  part  of  the  study.  It  should  be 
noted  that  the  above  array  has  90«  surface  wave  detection  capability  for 
ni|^  as  low  as  4.5  (Harley  and  Heiting,  1972).  However,  due  to  non- 
availability and/or  difficulty  of  getting  ALPA  data,  the  scope  of  this 
study  was  limited  to  a single  standard  station  (C.ollege)  data. 

DATA 

We  have  examined  College  station  long-period  data  from  1968  through 
1972  for  shallow  earthquakes  and  underground  nuclear  explosions  located 
in  central  Asia,  Nevada  and  central  Aleutian  areas.  Explosions  with 
mjj<5.7  are  not  recorded  by  the  long-period  components  (magnification 
M500)  near  T^20  sec  for  epicentral  distances  A>2000  km.  This  factor 
further  reduced  the  scope  of  the  study  and  we  have  based  the  investi- 
gation from  events  located  in  the  Nevada-Cal ifornia  and  Aleutian  regions. 

Initially  it  was  planned  to  use  central  Alaskan  short-period  array 
data  to  calculate  the  mean  body  wave  magnitude  m^^  from  multi-station 
values.  However,  the  stations  were  not  calibrated  as  frequently  as 
desired  due  to  limitations  of  funding.  Since  large  uncertainties  may  be 
introduced  in  nij^  computations  from  this  source,  we  considered  short- 
period  amplitude  data  only  from  the  College  station. 
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Evcrnden  (1967)  showed  from  United  States  short-period  data  that  the 
nijj  relation  depends  on  the  particular  phase  which  arrives  first  at  a 
given  epicentral  distance.  He  developed  magnitude  relations  for  a 
number  of  P-wave  phases.  As  discussed  in  Section  1 of  this  report,  the 
travel -time  studies  with  central  Alaskan  data  yielded  two  P wave  phases 
for  700  £ A <_  2800  kin.  The  apparent  velocities  obtained  for  these  tv;o 
phases  were  8.29  kni/sec  and  10.39  km/sec.  Since  appropriate  relations 
for  ni^  are  not  available  for  this  area,  we  considered  the  following 
two  expressions  given  by  Evernden  (1967) 

% = log(A/T)+2  logA-3.27  A<18°  (1) 

% = log(A/T)+4  logA-10.35  A>18°  (2) 

In  the  above  expressions,  (A/T)  is  the  ground  displacement  in  millimicron 
per  sec  and  A is  the  epicentral  distance  in  kilometers.  Equations  (1) 
and  (2)  were  proposed  by  Evernden  (1967)  for  P-wave  phases  which  ere 
associated  with  apparent  velocities  of  3.5  km/sec  and  10.5  km/sec, 
respectively.  These  two  P-wave  phases  seem  to  correspond  to  those  obtained 
from  central  Alaskan  data  for  700  ^ A < 2800  km  as  discussed  above. 

For  surface  wave  (T-v20  sec)  magnitude  computation,  we  have  used  the 
following  relations  (Nuttli,  1973) 

Mg  = log(A/T)+1.66  logA+2.6,  A<20°  (3) 

Mg  = log(A/T)+1.66  logA+3.3,  A>20°  (4) 

where  (A/T)  is  the  zero-peak  ground  displacement  in  milli  micron  per  sec 
and  A is  the  epicentral  distance  in  degrees. 

RESULTS 

Despite  very  careful  data  reduction,  the  body  wave  magnitude  (m^^) 
computation  by  Equations  (1)  and  (2)  resulted  in  values  which  showed 
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difforences  in  the  range  of  + 0.5  to  + 1 from  those  reported  in  NOAA's  Monthly 

Bulletins  (PDE).  It  should  he  noted  that  the  reported  values  given  in  PDE 

are  means  of  multi -station  measurements.  The  above  discrepancy  led  us 

to  conclude,  as  pointed  out  before  by  Liberman  and  Pameroy  (1969)  and 

otliers,  that  single-station  m^^  determination  has  inherent  large-scale  uncertainties. 

Thus,  we  used  NOAA's  mj^  values  to  investigate  the  relations. 

The  computation  of  surface  wave  magnitude  (M^)  by  Equations  (3)  and 
(4)  showed  a maximum  difference  of  _+  0.2  for  those  events  reported  in 
the  r>DF.  For  three  central  Aleutian  events,  the  comparison  is  shown  below. 


M^(NOAA) 

M^(Single 

Cannikin 

5.7 

5.7 

Cannikin  Collapse 

4.9 

5.02 

Mi  1 row 

5.0 

4.96 

Since  NOAA's  PDE  provides  M^  for  few  events,  and  as  the  single-station 
determination  showed  reasonable  agreement  with  those  obtained  from  the 
multi -station  measurements,  we  used  values  determined  from  College 
station  data  by  Equations  1 and  2.  The  results  of  plots  for  the 

Nevada-Cal ifornia  and  central  Aleutian  areas  for  earthquakes  arj  explosions 
are  shov/n  in  Figures  1 and  2.  In  these  plots  the  relation  derived 

for  the  w -model,  and  those  of  Gutonberg-Richter  (1956)  and  Thirlaway- 
Carpenter  (1966),  are  also  shown. 

In  the  central  Aleutians  the  three  undergound  nuclear  explosions  show  clear 
separation  from  the  shallow  earthquakes  in  the  same  area.  The  relative 
shift  on  the  left  of  the  Gutenberg-Richter  line  for  the  Cannikin  collapse  is  in 
close  agreement  with  those  reported  (Evernden  et  al.,  1971)  for  the  Benham 
explosion  and  collapse.  The  above  authors,  by  comparing  the  amplitude 
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spectra  for  tlie  explosion  and  collapse,  showed  that  the  collapses  are 
associated  with  a greater  proportion  of  long-period  energy  compared  to  the 
explosion.  In  tlie  relation,  this  implies  that  for  a given 

level,  the  corresponding  should  show  a relative  shift  toward  the  left 
and  is  in  agree'rient  with  the  results  in  Figure  1 and  2. 

The  Cal  ifornia-flevada  earthquakes  and  underground  nuclear  explosions  in 
Nevada  show  a clear  separation.  Moreover,  the  points  in  the  versus  m^ 
relation  for  this  area  in  comparison  with  those  for  the  central  Aleutians  seem 
to  indicate  a relative  sin  ft  to  the  left  for  a given  level.  The  cause  of 
this  feature  is  possibly  related  to  the  differential  path  attenuation 
(Ward  and  Toksoz,  1971);  this  needs  further  study  for  a complete 
explanation  of  the  above  phenomenon.  Hov/ever,  it  may  be  noted  that  in 
the  summary  of  Aki  (1972  , Fig.  7),  similar  results  as  above  v/ere  shown 
for  the  United  Stales  in  comparison  to  the  Aleulian-Kamchatka  area. 

CONCLUSIONS 

In  the  present  study  we  have  demonstrated  for  both  the  areas  considered 

a clear  separation  in  the  versus  relations  between  underground 

explosions  and  earthquakes.  We  have  also  considered  the  known  M -m. 

s b 

relations  (Gutenberg-Richter,  Thirlav/ay-Carpenter  and  Aki's  u^-model ) 
in  this  study.  So  far,  a best  fit  between  these  relations  and  the  present 
data  could  not  be  obtained. 
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riGURC  CAPTIONS 


Figure  1.  Ms-nij^  relations  for  underground  explosions  in  Nevada  and  earthquakes 
located  in  Nevada-Cal i fornia  regions. 

Figure  2.  Ks-m^  relations  for  undergound  explosion  in  Central  Aleutian  and 
earthquakes  from  the  neighboring  regions. 
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Figuro  2. 


Ms-m,  relations  for  underground  explosion  in  Central  Aleutian 
and  earthquakes  from  the  neighboring  regions. 
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